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ABSTRACT 

We use the perturbative QCD methods of Lepage and Brodsky to calculate 



> 
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r-fehe rate for pKsi with an eye toward the CP violating unitarity triangle 

m 

igle 7. We show that , although the penguins are large, there are regions 



00 

Qsf the allowed parameter space of the Cabibbo-Kobayashi-Maskawa mixing 
' 

mifitrix wherein 7 is measurable in the sense that penguins change the value 
I 

^ sin (27) one would extract from the attendant time dependent asymmetry 
^6asurement by less than 29%, so that a 3cr measurement of sin(27) as being 
different from is allowed by the corresponding theoretical uncertainty. This 
"would establish CP violation in Bs decays. The rates which we find tend to 
favour the type of luminosities now envisioned for hadron-based B-factories. 
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Now that there are two asymmetric e~^e~ colhding beam B-Factories, the 
SLAC-LBL-LLNL and KEK Asymmetric B-Factories, as well as several other 
B-Factory type machines, such as HERA-B, the CESR upgrade, and the Teva- 
tron upgrade, for example, under construction, the systematic exploitation of 
these machines for CP violation studies is not far away. To realize the true 
potential of these studies, it is important that the complete set of Standard 
Model CP violation parameters for the B-system be explored, if it is at all 
possible. In particular, this means that all CP violating angles a, (3 and 7 of 
the unitarity triangle should be measured ,where we use the notation of 
for these angles. The angle (3 is the "gold plated" angle of the triangle, as it 
will be presumably the most readily measurable of the three angles, via the 
modes B / JKs-,'^ / JK*^. It {13) is in fact used to specify the minimal 
requirements for the B-Factory machine and detector system to be successful. 
(Here, K*j^ denotes the CP + neutral K* meson.) Accordingly, the B decay 
modes needed for measurement of the angles a and 7 must also be identified 
and assessed. In this connection, the mode Bg^p + Ks is worthy of some 
attention; for, were it not for the possible contamination from penguins, this 
mode would be a candidate mode for the measurement of 7 Indeed, the 
potential contamination from penguins is just as substantial as it is for the 
mode B 7r%^ in connection with the measurement of a, for which the au- 
thors have devised isospin methods to combine the measurements of the 
modes B —>■ tt+tt^, tt'^tt" and B^ Tr+Tr*^ to extract a independent of the 
size of the penguin contamination- the main experimental problem of course 
is the measurement of the tt'^tt^ mode. It is desirable to address these pen- 
guin CP violation pollution effects from a dynamical approach which aims to 
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quantify them directly, thereby isolating just where a measurement may still 
be made, in view of the available parameter space in the respective Cabibbo- 
Kobayashi-Maskawa mixing matrix. Indeed, in a recent paper , we analysed 
the theoretical expectations for the size of these penguins in the basic mode 
B —>■ 7r+7r~ as well as in the companion mode B ti^ti^. We have found that, 
in a large region of the parameter space, the Asymmetric B-Factory devices at 
SLAG and KEK will be able to extract the fundamental CP violating angle a 
without depending on the penguin trapping methods in Ref. 0. The natural 
question to ask is whether an analogous region exists in the case of the mea- 
surement of the angle 7 in the Bg — > pKs decay? It is this question that we 
address in the following theoretical development. 

Specifically, we will use the approach of Lepage and Brodsky as it is 
represented in our analysis of ^ tt+tt^, K^K" in Ref. [Q. In this realiza- 
tion of perturbative QCD for hard exclusive processes, as we shall illustrate 
explicitly below, the exclusive amplitude is represented as a convolution of a 
hard scattering kernel (referred to as Th in Ref. [Q) with distribution ampli- 
tudes that sum the respective large QCD coUinear logarithms associated with 
radiation from the external legs of the constituent partons. These distribution 
amplitudes therefore obey a rigorous QCD evolution equation derived from 
QCD perturbation theory in Ref. ||^. We refer to this representation of hard 
exclusive hadron processes as the Lepage-Brodsky method. It was already 
formulated in Ref. in the context of the exclusive two-body B decays to 
light mesons of the type of interest to us here. See also Refs. for further 
illustrations of the method we shall use. As we explain in Ref. we expect 
the accuracy of our methods as used here to be at least as accurate as the 25% 
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accuracy determined in the work in Ref. We present both the absolute 
decay rates and the ratio of branching ratios corresponding to such rates, with 
and without the penguins included in the respective calculations. In this way, 
we expect to minimise the sensitivity of of our results to the uncertainty of the 
normalisation of the distribution amplitudes which we do use. Indeed, in the 
respective CP asymmetry parameter sin (27) analysis, we compute its appar- 
ent shift away from its expected value in the absence of penguins in ratio to 
that expected value, Asin(27)/ sin(27). We refer to this shift as the penguin 
shift of sin (27). The analog of this shift plus unity was already introduced 
in Ref. in the study of the time-dependent CP-violating asymmetry in the 
TT+TT^ decay mode. We will exhibit a formula for the penguin shift of sin(27) 
here for definiteness in complete analogy with what we have already published 
in Ref. for the corresponding shift of the analogous CP violating asymmetry 
parameter sin(2a) for the tt+tt^ decay mode. Evidently, the normalisation of 
our distribution amplitudes also drops out of the penguin shift of sin(27). 

Concerning the Cabibbo-Kobayashi-Maskawa (CKM) matrix itself, we fol- 
low the conventions of Oilman and Kleinknecht in Ref. [jlO| for the CP- violating 



phase (5i3 = 6 and in view of the current limits on it we consider the entire 
range < 5 < 27r. For the CKM matrix parameters Vtd and Vub we also 
consider their extremal values from Ref. |T0| (the Particle Data Group (PDG) 
compilation). To parametrise these extremes, we use the notation defined in 
Ref. in|] for \Vub/Vcb\ in terms of the parameter Rb = .385 ± .166 0]. All 
other CKM matrix element parameters are taken at their central values 



We should emphasise that the decay under study here is not the only way 
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to study the CP-violating angle 7. Indeed, due to the very small rates which 
we shall find, it will be seen that the most appropriate machine to pursue the 
mode under discussion here is a hadron collider type B-factory device. As 
shown in Refs. |TI], |12| , the e^e~ colliding beam type B-factory device can 



approach 7 from other decay mode avenues. 

We further emphasise that it is possible to use the methods of Lepage and 
Brodsky ,as they are represented in the analyses in Refs. P,^, to address 
both the concept of colour suppression for the —* pKs decay as well as the 
size of the penguin pollution in its CP violating phase structure as described 
above. We will take advantage of this opportunity to get a quantitative es- 
timate of the colour suppression effect in this decay under study here. In 
practice, what this will mean is that, in addition to computing our branching 
ratio (BR) for the decay with and without penguins included, we will also 
compute it with and without gluon exchange between the would-be spectator 
s and the qq lines of the outgoing p. Again, we will focus on the respective 
ratios of BR's to avoid sensitivity to the uncertainty in the normalisation of 
our distribution amplitudes. Such an estimate of colour suppression has not 
appeared elsewhere. 

Specifically, we note that the QCD corrections to the weak interaction La- 
grangian will be represented via the QCD corrected effective weak interaction 
Hamiltonian ?^e// as it is defined in Ref. |p | 



{2 10 
fc=l k=3 



+ h.C. (1) 



V2 

where the Wilson coefficients Ci and operators Qk are as given in Ref. |Tl 



Gp is Fermi's constant, p is is the renormalization scale and is of Oirrib) and 




Figure 1: The process Bg ^ p + Ks- The four-momenta are indicated in 
the standard manner: Pa is the four-momentum of A for all A.To leading 
order in the perturbative QCD expansion defined by Lepage and Brodsky 
in Ref. 0, the two graphs shown are the only ones that contribute in the 
dominant contribution as isolated by the methods of Ref. when penguins 
and colour exchange between the outgoing p partons and the outgoing Ks 
partons are ignored. The remaining graphs in which the gluon G is exchanged 
between the would-be spectator s and the remaining p parton lines as well as 
the penguin type graphs are shown in Figs. 2 and 3, where we see that, for 
QCD penguins, there is the added possibility that the gluon G interacts with 
the penguin gluon itself of course. 

here q = s. The application of this effective weak interaction Hamiltonian 
to our process Bg pKs then proceeds according to the realization of the 
Lepage-Brodsky expansion as described in Ref. @]. This leads to the "domi- 
nant" contribution in which the p is interpolated into the operator O2 = Qi 
in Heff via the factorised current matrix element < p|'u(0)7^Plm(0)|0 >, 
Pl = |(1 — 75) so that the respective remaining current in O2 = Qi is re- 
sponsible for the Bg to Ks transition shown in Fig.l. We refer to this con- 
tribution as the "Tree" contribution. The complete amplitude for the process 
under study here, Bg pKs, is given by the sum of the graphs in Fig. ^ 
and those in Figs. ^ and ^, to leading order in the Lepage-Brodsky expan- 
sion defined in Ref. and realized according to the prescription in Ref. 0. 
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Figure 2: The colour exchange graphs for the process Bs p + Ks to leading 
order in the Lepage-Brodsky expansion in Ref. 0,0, ignoring penguins. The 
kinematics is as defined in Fig. 



In Fig. ^ we show the graphs in which colour is exchanged between the 
would-be spectator s in Fig. |1| and the outgoing p parton lines and in Fig. ^ 
we show the respective penguin graphs: the dominant graphs according to 
the prescription in Ref. (3a, 3b), the colour exchange graphs (3c, 3d), and 
the exchange of the hard gluon G between the would-be spectator s and the 
penguin gluon itself for QCD penguins, (3e), which we also will classify as 
colour exchange. To address the issue of factorisation/colour-suppression, we 
shall present results when graphs in Figs. 2 and 3c-3e are dropped and when 
they are included. We thus give results for the approximations in which only 
the graphs in Fig. |l] are included (Tree), in which the graphs in Figs. 1, 3a 
and 3b are included (Tree+Penguin), in which the graphs in Figs. 1,2,3a and 
3b are included (Tree+Penguin+Tree Colour Exchange(Ci?r)), and in which 
all graphs in Figs. l,2,3a-3e are included (Tree-|-Penguin-|-Tree and Penguin 
Colour Exchange(C£'r+p)). For the electroweak (EW) penguins, there is no 
penguin gluon with which the would-be spectator s could interact. We need 
to stress that, as shown in Ref. the usual QCD factorisation properties for 
exclusive amplitudes at large momentum transfer are sufficient to justify the 
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(e) s ^Ks 

Figure 3: The penguin graphs for the process Bg —>■ p + Ks, to leading order 
in the Lepage-Brodsky expansion defined in Ref. 0,0. The kinematics is as 
defined in Fig. |^. 
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formulation of our amplitude according to the graphs in Figs. 1-3. More phe- 
nomenological arguments, such as the current field identity based BSW model 



in Ref. |jT3[, etc., which would lead to the same graphs, are not needed. 



Some discussion of the effective values of the coefficients Ci = C*2, 6*2 = Ci 



in relation to the coefficients ai and 02 as defined in Ref. [13| is now appropriate 



Following Ref. |]T3[ and the recent results in Ref. |jT^, when we use the standard 
QCD to calculate the diagrams in Fig. 1 and take them alone as our result ( 
this is our definition of factorisation) , we use 02 = .24 = |C2(mfe)| and when 
we assess the colour-suppression effect by including the exchange of G between 
the s and the qq of the p we set Ci(mfe) = 1.1; these results are consistent with 



those found in Ref. [Q. We note that the naive relation 02 = C2 + |Ci = .127 



would give a value for 02 that is about a factor of two smaller than what 
is found in Refs. [|T^, and the references therein. The parameters ai,a2 
are therefore purely phenomenological properties of the hard effective weak 
interaction process and can be taken from experiment in our analysis: one 
may view 02, for example, as the effective value of C2 + \Ci when the current 
field identity is used to interpolate the p into our effective weak interaction 
vertex. The Lepage-Brodsky formalism then allows us to calculate the recoil 
corrections associated with the momentum transfer required for the would-be 
spectator to be kicked from the Bg to the final outgoing Ks using perturbative 
QCD to describe the respective hard gluon exchange, as we noted above. This 
"kick" is the defining aspect of our calculation of Bs —>■ pKs in comparison to 
those in Refs. |rS] and in fact in comparison to the related two body B decay 



analyses in Refs. [|T^]. The point is the following. As one can see from the 
results in Refs. contrary to what happens in the tree level part of the 
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calculations in Refs. [0,|TB], the graph in Fig. la in which the hard gluon kick 



to the spectator comes from the b-quark line (the heavy quark line) develops 
an imaginary part that is treated rigorously in our work so that there is a 
non-trivial strong phase for our "tree level" contribution compared to those 



in Refs. [15,16]. This happens because, as > rub + rus where niq are 
evaluated at the scale ^ rriB, the heavy quark line can reach its perturbative 
QCD mass shell in the graph in Fig. la, and in the similar graphs in Figs. 2 and 



3. Evidently, this effect is missing in the results in Refs. |]T5|,|T6|. Any serious 



discussion of the CP asymmetries in the amplitudes for exclusive two-body B 
decays must take this strong phase into account in general (it is different for 
Tree and Penguin contributions for example) as one can see from our formula 
for time- dependent asymmetry in Bs — > pKs below. Our paper is the first 
paper to do this systematically. 



Here, we should also comment on the recent results of Ref. 17 on the 



process B ^ nn. The authors in Ref. jl^ use the same Feynman diagrams, 



analogous to those in Figs. 1-3 here, as we have shown already in Ref. 
and same Lepage-Brodsky expansion formalism except that they assume the 
graphs analogous to those in Fig. 1 are to be replaced by a real form factor 
with the appropriate external wave function/decay constant factors. The usual 
corrections to the diagrams in Fig. 1 are then represented as a power series 
in as times this assumed real form factor. We do not make such an assump- 
tion; we calculate systematically in the Lepage-Brodsky expansion. A major 
difference is that the authors in Ref. miss the recoil phase of the dominant 



contribution in the analog of Fig. 1 for the B ^ nn process, although they 
do calculate the recoil phase in the respective analogs of Figs. 2 and 3. To 
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see what effect this has, we note that, from our Eq.(5) in Ref. we get the 



direct CP violation result for the B ^ nir process as 



-0.0086 < A^Sp < 0, for 7 G (0, vr) 

-0.0086 < Acp < -0.0050, for ^ > 7 > ^, (2) 

with Aqp monotone decreasing in the second currently preferred region of 

7 for 7 < 7 < 1.806 and monotone increasing for 1.806 < 7 < whereas in 



Ref. [|T^ this asymmetry is predicted to be —4% x sin 7. Evidently, experiment 



will soon be able to distinguish between the two approaches. See Ref. [IS] for 
further discussion of this and related matters. 

In this way, using the methods of Ref. [Q we evaluate the graphs illus- 
trated in Fig. 1-3 and arrive at the results in Table 1 and in Fig. 4 ( the 
explicit expressions for the respective amplitudes may be inferred from those 
for the process B ^ mi given in Eq.(l) and in Eqs.(Al-A4) in Ref. [§] via the 
appropriate substitutions of momenta and distribution amplitudes; for exam- 
ple, for the factor Fj^ in (Al)in Ref. we would now have its form obtained 
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by the substitutions 



-^yw2 -^^ymii + 3/? Km - yi)) 

P-K+ Prs 
rriu ma 
rriKs 

Axly2El+ml q^P^ + qf - + le] 
{-2x2y2E^+mB){yim\ - ml + ie) q\Pb - qf - ml + ie] 

(3) 

wherein q = Pg — Pg, Pf, f = b, s is 4- momentum of / in the Bg in Fig. 1 
and PdjPg are the 4-momentum of d,s respectively in the Ks in Fig. 1, so 
that we have Ps = X2Pb^ and P^ = y2PKs^ for example, and [3' ^ — -418 is the 
asymmetry parameter in the Lepage-Brodsky distribution amplitude for the 
Ks as determined in Ref. PD[ and evolved to the scale m^. We use fx — 0.112. 



In this regard, we further note that the Lepage-Brodsky distribution amplitude 
for the p in the analog of Eq.(A4) in Ref. for the process under study here 
would substitute VSfpTUp /{Pp)ziZ2{lM8 - 1.7Azi + 1.7Azf) for 

fTrZiZ2'j5{ f^o + rriT,) for example by the standard methods, where we use the 



Chernyak-Zhitnitsky (C-Z) type result |21] for the p distribution amplitude 
in analogy with our discussion in the Notes Added in Ref. 0. Here, e(Pp) 
and fp are the respective p polarisation 4-vector and decay constant with 
fp = .lAGeV. The Bg distribution amplitude is taken in complete analogy 



with the Bd in Ref. 0, so that it is given by aB(f)B{'Wi, W2)/ \/2Nl = aB^{w2 
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X2)l^/2Nc where Nc = is the number of colours, = fsJ^/^Nl and X2 = 
0.0542 is determined, as we present in our Appendix, following the treatment 
of heavy mesons in Ref. 0] using potential model parameters such as those 
in Ref. ||22|. Finally, note that the quark masses are the running current 
quark masses For completeness, the complete result for the amplitude 

corresponding to the graphs shown in Figs. 1-3 is given in the Appendix. 
Moreover, the precise definition of the penguin shift A sin(27) is given by the 
following generalisation to our process Bg pKs of the formula of Gronau in 
Ref. for the corresponding shift of sin(2a) due to penguins in the — » vrvr 
process 

- sin(27) - A(sin(27)) = ^ (4) 

for 

)^ ^ ('51 

ATe+''t'T+^ST + ^ . Ap^e^'^^'J^'^'-^ ' 

where the amplitude Atc'^'^'^^^^'^ corresponds to the tree-level weak processes 
in Figs. 1 and 2 and the amplitudes Ap^e~'''^^^~^^^^^ correspond to the respective 
penguin processes in Fig. 3. Here, we identify the weak phases of the respective 
amphtudes as 0.^, r = T, Pj and the attendant strong phases as 5^, r = T, Pj. 
In general, j = 1,2 distinguishes the electric and magnetic penguins when this 
is required, as one can see in our Appendix. In this notation, we have 7 = 0t- 



From the results in Table 1, and their ratios with one another, we see that the 
colour suppression idea does not really hold for this decay. We see, as already 
anticipated by several authors [jl|, that the penguins are indeed important. 
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1.000 








1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


L 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


J 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




0.500 














0.000 






















-0.500 










-1.000 


1 1 1 1 li 1 


1 1 1 1 1 1 1 


ill Ill 




1.000 


2.000 3.000 4.000 


5.000 6.000 



Figure 4: Penguin shift of the CP asymmetry sin(27) in Bg — > pKs for 
Rj, — 0.385 for the matrix element approximation corresponding to the last 
column in Table 1. The analogous plots for the ilcr values of Ri, are discussed 
in the text. 
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[0.00624,0.151] 
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Table 1: BR for Bg — > pKs as a function of Rb as defined in the text. The 
factorised approximation without penguin effects is denoted as Tree; the corre- 
sponding results with the penguin effects (both EW and QCD penguins) included 
are denoted by Tree + Penguin; the results corresponding to the inclusion of 
the gluon exchange between the uu in the p and the s would-be spectator are 
denoted by Tree -\- Penguin -\-CEt; and, when the gluon exchanges between the 
s would-be spectator and the outgoing dd of the p and the penguin gluon itself 
are included, we denote the result by Tree + Penguin + CEt+p- All results 
are given with a factor of (/B^/.141G'eV^)^ x 10"^ removed for a total width 
T{Bs all) = 4.085 x IQ-^^GeV and for the variation 0.0 < 5i3 < 2tt. 

There is a regime, 0.0" < 7 < 40.5°, 102.5" < 7 < 157.9°, for the central 
value of Rfy for example, wherein the shift of sin(27) is less than 29% of its 
magnitude so that it would be measurable in this regime if the luminosity 
is large enough to provide a sufficient number of events. By measurable, we 
mean that a 3a result for its value is not blocked by the uncertainty from 
penguins. We define this regime in which |Asin(7)/ sin(7)| is less than 0.29 as 
the measurability regime. Approximately 34% of this regime of measurability 
intersects the allowed region given by the limits on 7 discussed in Ref. [jl|, 
135° ^ 7 ^ 45°. We need to stress the following. When the pollution in 
the sin(27) is < |sin(27)|, a 15 — 20% accuracy calculation of the pollution 
is sufficient - when we have as well |A(sin(27))| < .29|sin(27)| sin(27) is 
directly measurable; when .29| sin(27)| < |A(sin(27))| < | sin(27)|, we measure 
a quantity from which we can extract sin(27) with 20% theoretical precision 
so that sin(27) can still be extracted. However, when the pollution is itself 
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dominant and sin(27) is ~ 20% of it, a 20% accuracy knowledge of the pollution 
will not permit the extraction of sin(27). Thus, for a given precision on the 
pollution, depending on the relative size of the pollution and sin(27), one has 
these three regions and one of these is exactly that addressed as our regime 
of measurability, one wherein sin(27) is measurable. (For the ±lcr deviations 
of Rb, the measurability regimes are qualitatively similar in size and location, 
with the exception that the lower regime is absent for the —la deviation case. 
So, we do not show these ±la deviation measurability regimes separately here- 



see Ref. |]T9| for the corresponding plots analogous to that in Fig. 4.) 



The question naturally arises as to the sensitivity of our regime of measur- 
ability to the parameters in our calculation. We now turn to this. Since the 
penguin shifts plotted in Fig. 4 are determined from amplitude ratios, they do 
not depend on the normalisations of the distribution amplitudes, or the hard 
recoil gluon exchange coupling in Figs. 1-3. What they do depend on are the 
relative strengths of the leading and non-leading Gegenbauer coefficients @] 
in the distribution amplitudes, the relative strengths of the penguin and non- 
penguin operators in the effective weak Hamiltonian ( 02 and the value of 
in our one-loop penguins), the quark running masses and the light-cone frac- 
tion X2 of the s in the Bg as determined by the Cornell model of B mesons. 
We have varied all of these parameters systematically as currently allowed by 
the la limits on them when they are taken from data or theory together with 
data [|19| . We find that the ffist part of the regime of measurability varies from 
[0°, 19°] to [0", 58"]. Thus, it may even be true that some of the allowed regime, 
45" ^ 7 ^ 135, overlaps this ffist part of our regime of measurability. The 
most important aspect of this variation is that most of it is due to changing the 

15 



value of the running b-quark mass by just ±3.5% and by varying the value of 
02 between .14 and .34 (for reference, the variation in as is just that generated 
by the la variation of Aqcd (see the following), the variations of the non- 
leading Gegenbauer coefficients are the la variations as determined from their 
extraction from data in Refs. [2^,21], the la variations of the running quark 



masses are as given in Ref. |^ and the methods therein , and the variation of 
X2, between 0.041 and 0.071, is as given by the parameter variations allowed 
in Refs. - see Ref. for further details). If we do not vary these two 
measurable parameters, then the first part of our measurability regime only 
varies between [0", 35°] to [0°,47°], i.e., it is robust to the remaining parame- 
ters in our calculation. In the actual precision hadron B-factory environment, 
we can expect that both rUb and 02 will be known much better than we know 
them currently from comparison with data, either experimental or theoretical 
(lattice) data. The current large sensitivity to and 02 of the upper bound- 
ary on the first part of our regime of measurability is mainly academic because 
this regime is already outside the preferred region of 7 and the variations we 
see with m^, 02 still leave most of this first part outside the preferred region. 
The second part of our regime of measurability begins at 102.5° and ends at 
157.9°. Upon variation of our fundamental parameters as we described above, 
the beginning point varies between [98.2°, 105.5°] and the ending point varies 
between [138°, 180°] , so that the preferred region of 7 which overlaps the sec- 
ond and most important part of our regime of measurability, [102.5°, 135°] is 
only changed by 3 degrees by the current uncertainties in our fundamental 
parameters. Again, if we do not vary 02 and rrih, this already small effect is 
reduced significantly. We thus have a robust prediction that 7 is measurable 
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in the regime [102.5^135°]. 



Recently, several authors have argued that current data actually prefer 
the regime 36° < 7 < 97°, although more recent theoretical analyses [^, ^ 
would question this conclusion. Here, we stress that, from our results in Fig. 
4, we can see that, in this new so-called preferred region, except for the small 
region 86.6° < 7 < 92.7°, the penguin shift is bounded in magnitude by a 
factor of 2 relative to the actual value of sin(27) so that, as we have a ~ 15% 
accurate knowledge of this shift, we still may use our results in the Appendix 
to radiatively correct this pollution out of sin(27) to the ~ 30% accurate level, 
allowing again a 3cr measurement of sin(27). The use of this technique to make 



fundamental tests of the SM is well-known p7|. 

The BR's in Table 1, which remain qualitatively similar to their values 
shown here under the variations of parameters just considered, however tend 
to indicate that the required luminosity would be more appropriate to hadron 
machines than to an e+e^ annihilation B-factory. We note that the results 
in Table 1 are somewhat lower than the general range of similar results in 
Refs. [|^,|l^]. For example, our highest values for the BR just reach the lowest 



values in latter references. The recent and upcoming measurements of rare B 
processes can then already discriminate among various models of these pro- 
cesses on the basis of decay rates alone. To illustrate this,we note that in 
Ref. we used our methods to compute the range 

1.87 X 10-^(^?,2(m|)/<7,2(m|)|^(.)^^„^^^^)2(/5V0.136GeV)^ 
<BR(5rf^7r+7r-) < (6) 

2.63 X 10-6((7f(m|)/(7,'(m|)|^(5)^^^^^^^)2(/5^/0.136GeV)l 
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(5) 

We note that, according to Ref. |£8[, the current two-loop value of ^qcd 
2?)7^^eV and according to Ref. the best value of V^/b^ is now 210 ± 
30MeV so that we have the estimate 

(^7s("^l)/^?'(^l)lAg)^=o.iGev)'(/^/0-136GeV)2 - 1.70. (7) 



This means that our result in Eq.(^ is consistent with the recent CLEO re- 
sult 13 BR{Ba ^ TT+TT-) = 4.7^^1 ± 0.6 x 10"^. Nonetheless, even if we 



allow the entire range which we and the authors in Refs. [|1^, [T^ find for 
BR{Bs pKs), we are led to suggest that the B-factory of the SLAC- 
LBL/KEK type should focus its attention on other possible roads to 7. Oth- 
ers |]1[ have reached a similar conclusion. 

Finally, we stress that we have found that the assumption of colour sup- 
pression (factorisation) does not appear to work very well in our calculations. 



This is consistent with the results in Refs. [^,|32| on the analysis of the data 
on the processes B ^ ^! / J K^*\ We will take up the corresponding analysis 
with our methods elsewhere [|l^ . 
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Appendix 



In this Appendix, we record for completeness the amphtude which we have 
evaluated from Figs. 1-3. Specifically, following the usual Feynman rules and 
the prescription given in Ref. |0, as already illustrated in Ref. 0, we get the 
amphtude 



2mB2Ep2EA's(27r)9/2 



(A.l) 

where the "would-be tree level" contribution to the amplitude is, from Figs. 1 
and 2, given by 



ciB(l)Blb{pB - rriB) 



V2V2 



V2 



I ■ ^c ^ fK(l)Klh{fKs+^Ks) , ■ .c c 



V2V2 



'-iGFa2VubVl 



-zGirasKfeK^^e^/,^^ _ ^^^■^ aB(pBl5{fB-mB) ^^ 



V2 

+ Tee / d[z\ 



v/2 



+ 



A'^TmII - 75) 



/p0p^(PpK,e^^^,_^^)) 



^2^2 



J r ) 



(A.2) 



where contribution of Fig. 2 is (not) included for rce = 1(0) and where 



^(f>K{yi,y2), ^(pp = 7373*^^(^1' ^2) and = -^^(j)B{wi,W2) are 
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the Lepage-Brodsky distribution amplitudes — 

Ml, 1/2) = ym (1 + 3/5^(1/2 - yi)), M^u ^2) = ^12:2 (1.348 - 1.742:1 + 1.74^2), 

and (f)B{uii,W2) = 6{w2 — X2) are as indicated in the text above with X2 = 
(m^ — [ml + ml — mB)ml/ (m^ + ml)) following the treatment of heavy mesons 
suggested by Ref. [ffl based on non-relativistic potential model considera- 



tions for example. Here, the constituent quark masses are taken as [22 



m'l = 0.51GeV and mf = 5.1GeV, so that X2 = 0.0542 when we take 



5.369GeV, as we should according to Ref. |]I0|. From Ref. we have 05 = 
/b/"\/12 where /b is the B decay constant. Here, Pa is the 4-momentum 
of A for all A and, when a parton-type occurs in two external wave func- 
tions a prime is used to distinguish the two 4-momenta in an obvious way. 
To be precise, let us list these internal parton momenta as follows for Fig. 
1 : P+ = XimB, P,- = {ml + Ql{B))/ (ximB), Pb± = Q±{B), P-^ = 
X2mB, Pf = {ml + Ql{B))/{x2mB), Ps± = -Q±{B), Ql{B) = xiX2m\ - 
X2ml-xml Pj = yiiEK+PKz), Pd = {m'i+Ql{K))/{y^{EK+PKz)), Pd± = 
Q^{K), P;+ = y2{EK + Pk.),P',- = K + Ql{K))/{y2{EK + Pk.)), PL = 
—Q±{K), Q\{K) = yiy2m\- — y2fnl — yim^, where we always work to leading 
order in Q^/m^ for all Qx, and where we use the usual light-cone notation 
with Ek = P^ and and Pkz = Pr so that P^ = P% ± P|., etc. The 
are the QCD colour matrices generating the vector representation carried by 
the quarks so that Qs is the QCD coupling constant. Thus, Eq.(A.2) illus- 
trates explicitly how the Feynman diagrams in Figs. 1-3 are evaluated for 
readers unfamiliar with the methods we used in Ref. P], for example. The 
standard trace and integration manipulations, taking into account the defini- 
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tion [dx] = dxidx25{l — xi — X2), then lead from Eq.(A.2) to the result 

Are-'^-e^'- = ^ T {V:,Vu,/CF){a,/V2){PcMsml/Q^) 
[hi {rriB - 2{mK + rUb) + mKmb/rriB) 
+ I22 {2imK - Ek) + m]^/mB) 

+ {x2mB - 2{x2Ek + XiUIk) +rnd + mKirriK - md)/mB) (0.291/L)2) 

+ (-4CiV^r,e/(mBa2))((0.291(-m| + ml 

+ mj)/(2m|))(1.348ai - 3.088a2 + 3.48a3 - 1.74a4) 

- (0.166666mJ/m|) (-0.25360 + 2.7586i - 2.50562))] 

(A.3) 

where the various mathematical symbols are defined below. Continuing in this 
way, using the entirely similar methods, we find that the penguin graphs in 
Fig. 3 correspond to the contributions to the amplitude in Eq.(A.l) given by 

Ap,e-'^^^e''^^ = ^ {{mBPcMs/Q^)ap'-''^himl[{x2 - 2{x2Ek + XimK)/mB + m\lm\) / D2 

V 2 

+ (1 - 2{mb + mK)/mB + m^mb/m|) /21/O.29I + (-1 + 2mK/mB + mJ/m|)/22/0.291] 

Ap.e-^'t'^^e''^^ = ^ {mBPcMs/Q^W\{mbml{x2mB/2 - mK)/D2) I33 

+ (^brriKXiiml - m\) + mbmB{{l + "^Xi - xl)m\ - X2m\ + (1 - ^Xi)x2m^)^ /32/-D2 

mbnipXi - -xim^m^ - 4xim6m^Ex(l - 7, ) + ^Ximbm^mKi- ) /32/21/O.29I 

2 2mB 2 rriB J 



+ ( 3ximbmBm% - XimbmK{m]^ + -m\ -m^) ) /32/22/O.29I 



1 o 1 

2 



- mbm^mBil - mb/ {2m b)) I 33 hi/ 0.291 
+ {mbmKml/0.582)h3h2]. 

(A.4) 
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In Eqs.(A.2)-(A.4), the following definitions have been used: 
:F^fKfpaBCj.glV3 

= {{Ek + PcMs)/{'^niB)){x^ml - ml + m^) 
721 = (-0.253/m|)4i + (2.505/m|)42 
I22 = (-0.253/m|)£22 + (2.505/m^)£23, for 
£21 = 0.403041 - 2.202003i 
^2 = -0.3794583 - 0.6585764i 

= -0.5097241 - 0.1969674i, 

731 = 0.0485 

732 = (0.291/.3)(0.195517/m| - 0.064303i/m^) 

733 = (0.291/.3)(0.132055/m| - 0.0608173i/m|) 

D2=ml + ml- {Ep + PcMs)ximB - {Ep - PcMs){x2mB + ml/ tub) 
ai = —1.0 — X2ln{xi/ X2) — X2m 
0,2 — —0.5 — X2 — a;2 ln(a;i/a;2) — xlTii 

03 = — 1/3 — X2/2 — x1-\- x\ \tl{xi/x2) — x^ni 

04 = -0.25(a;^ - xj) - (4/3) (x? + xl)x2 

— 3xl{xl — xl) — Ax\ — X2 ln(xi/x2) — X2ni 
bo = -ln(a;i/a;2) - ni 
bi = 61 

62 = 0,2 




22 



+ V;*\/i5{-4 (3a;? ln(x,)/(4(x( - l)^) + xt/A - 3xi/(4(a;, - 1))) 

- 2 (5V(2(a;t - 1))) (1 - Hxt)/{,Xt - 1)) - 2\VtdW-^xl\n{xt)/ {2{xt - if) 

- Xt{Q.2b - 9/(4(x4 - 1)) - ?,/{2{xt - 1)'))} + (4sy3)(0.641 - a;t(7/(3(x-t - 1)) 
+ U/{l2{xt - If) - l/{2{xt - If)) - xM^t){ll{^{xt - 1)) - 35/(12(xt - 1)^) 

- 5/(6(x, - If) + l/{2{x, - 1)')) + (2/3)2 ln(a; J) - (a;,/2 - 3/(4(a;, - 1)) 
+ 3(2x? - xt) Hxt)/{^{xt - If) - 0.75)} ] 
Xj = myM^, j = u, c, t, 
o^p^'^ = ^TO.(-0.195) + ^V;,V,, (0.641 + x,{l/(2(x, - 1)) 
+ 9/(4(x, - 1)2) + 3/(2(x, - 1)3)} - 3x3ln(x,)/(2(x, - 1)^)] 
Vcea = {CG/CF){mBPcMs/{4:Xc2dbkdbp)) [^cpoo{(-^ + 4a;i)m|- + (6 - 4xi)ml 



- 3xim| - llximftrmB + ^("^^s ~ "^d)} + ^cpio{5m| - Sm^ - GmBrriK} 

5 3 3 

+ ^{2^4 - + 2"^^^ (2m|/4^)[a;iidpoo{(a;imB + mK)EK - (2 - a;2)m^ 

777.2 12 1 

- Pp- Pk- X2mBmK + ^mxEp H ^(1 - 2mK/mB) - m^(- H (1 - 2mK/mB)) + ^^(m^ + 

3 777- 

+ idpoi{2EpEK + 2Pp ■ Pk{xi - mK/{2mB)) + mKEp{l - 5xi - 2~) + '^xitubEp - ^^'^s 

- 2ml{x2 + - ^) + (^ + A)m2} - a;iz,,io{2mB^x(l - ^) + m^} 



2ximQ niB Xiiri^ 2ms 2mi 
— 2 TidpnPp • Pk] 

XiUIq 
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2,2 2 

— 171^ + rrij^ — nip 


J 


2 2 2 

= rriB - rriK - rUp 


J 

dbr 


2,2 2 

— -\- nip — nij^ 


Xc2 


— /dbk 


^cpOO 


— u.ozy 






icpOl 


— Oa/o + 06/4 + Oc/o, wnere 


f0 


= U.41(5 


r 

Oa 


1 O /I o 

= l.o48 


Sb 


= -1.74 


Sc 


= 1.74 
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idpOl — ^airdp2{zo) + Sbirdpsizo) + Scirdp4:{zo) 

Wo = i^airdpiizo) + hirdp2{zo) + 5cVdp3 (^o)) (1 " rf})/2 = idpoo(l - r/3)/2 

idp02 — ^airdp3{zo) + ^birdpiizo) + ^cirdpbizo) 

idpll = {Sairdp2{zo) + Sbirdp3{zo) + Scirdp4{zo)){l - rp)/2 = idpOl(l - rp)/2 

idpi2 = (Sairdpsizo) + Sbirdp4{zo) + Scirdp5{zo)){i " rp)/2 = idp02{l " rp)/2, where 
irdpiiz) = -1 - 2;ln((l - z)/z) - nzi 

irdp2{z) = —.5 — Z — z'^ — z)/z) — TTz'^i 

irdpsiz) = -((1 - zf + z^)/3 - 3z{{l - zf - z^)/2 - 2,z^ - ^Mn((l - z)/z) - nzH 
trdp4{z) = -((1 - z)' - z')/A - Az{{l - zf + ^3)/3 - 2,z\{l - zf - z^) 

-^z^ - z^hi{{l- z)/z)-t:zH 
^rdp,{z) = -((1 - zf + - 5z{{l - zf - z')/A - lQz\{l - zf + z^)/^ 

- hz^{{l - zf - z^) - 5z^ - znn((l - z)/z) - ttzH, for 
Zq = ml/{xidbr) 

Vceb = ((1 - .5CG/CF)(mBPcMs)/(2Q^))(-l + rf3){Sairdpl{xi) + {Sb - Sa)irdp2{xi) 
+ {Sc - Sb)irdp3{Xi) - Scirdp^iXi)}, 

(A.5) 

where the kinematics is the usual two-body decay one: uib = + E^, E^ = 
dbkl {2mB), PcMS = ^A(m|,m^,m2)/ {2mB) for A(a;, y, z) ^ + + z'^ - 
2xy — 2xz — 2yz, so that the decay width itself is given by 

r{B ^ pKs) = \M\''PcMs/{S7rml) 

(A.6) 
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Here, Cg = 3, Cp = 4/3 and we have used the values |^ m„(lGeV) = 
S.OMeV, md(lGeV) ^ 8.9MeV, m,(lGeV) ^ .175GeV, m^{m^) = 1.3GeV, 
nibimb) = 4.5GeV, and mt{mt) = 176GeV. We take = sin^^iy = 0.2315, 
where 9w is the usual weak mixing angle; and, Oem is the QED fine structure 
constsnt. We note further that we use an average value of the square of the 
momentum transfer to the would-be spectator in Figs. 1-3 to get = 3.72 
in JF above; the analogous average for the square of the momentum transfer 
through the penguin yields — -25 in the evaluation of coefficients a^p above. 
Thus, in both cases, we see that the momentum transfers are large enough that 
they are well into the perturbative regime where the methods of Ref. [^] apply. 
This completes our appendix. 
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